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We study possible contributions to the lvtt'^ and Df — > p^vr^ decay amplitudes. The Df ujn'^ 

decay amplitude vanishes when naive factorization is used, while the — » p^vr"*" decay amplitude arises due to 
the annihilation contribution. We find that amplitudes for both decays might be a result of the internal K, K* 
exchange. The Df lutt'^ amplitude might obtain additional contributions from Df P'^viv') re-scattering. 
The low experimental bound on the Df — > p^Tr"*" rate can be understood as a result of combination of the 7r(1300) 
pole dominated annihilation contribution and the K, K* internal exchanges. The calculated branching fractions 
for Dt (j-'Tv'^ and Df — > p^Tr^ are in agreement with the current experimental results. 



The weak nonleptonic decays of charm mesons 
were usually approached within the factorization 
ansatz |ll2i;-{l4lf^ . A decade ago it was real- 
ized |2I3| that one has to include the effects of 
final state interactions (FSI), with the simplest 
approach being to treat the FSI by assuming the 
dominance of nearby resonances. This leads to 
rather good overall agreement with the experi- 
mental data 2 3 ; however, there are a few cases 
where none of the existing approaches work. Two 
such examples are the channels (quoting the PDG 
experimental values [2]) 

BR{Df -> ujn+) = (2.8 ± 1.1) x IQ-^, 
BR.{Dt ^ p\+) <7 xlO-\ 

The current theoretical approaches usually pre- 
dict that the Df p°7r+ branching fraction is 
equal [3] or even larger than the branching frac- 
tion for the D'^ wtt"*" decay 5 6^ in contradic- 
tion with the present data (Q). 

On the other hand, the observation of the 
Df W7r+ decay ^ has been motivated as a 
clean signature of the annihilation decay of 

The sizes of annihilation contributions are 
very important for phenomenological studies, but 
are also very hard to obtain from theoretical con- 
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siderations (see e.g., Pj). Understanding the ori- 
gin of the Df WTT"*" transition is thus of great 
theoretical interest. 

Let us first discuss the two modes ^ using fac- 
torization approximation for the weak vertex. In 
this approximation the — > ton'^ amplitude is 
zero due to G - parity conservation, which gives a 
vanishing {LUTr^\{ud)v-A\0) matrix element |1(J) . 
The p°7r+ decay amplitude, on the other 

hand, already in the factorization limit receives a 
contribution through the annihilation graph. Fig. 
1, 

M{Dt - p\+) = ^K.K:,aiX 

x(p\+\iud)v-A\0)msc)v-A\Dt), 

leading to simple tt pole dominance in the 
{p^Tr^\{ud)v -a\Q) matrix element. The analysis 
of [E] indicates that 7r(1300) states dominate this 
annihilation graph, while the contribution of the 
lowest lying tt is negligible. In [ni] we have es- 
timated the size of the annihilation contribution 
coming from the 7r(1300) intermediate state. We 
found /Tr(i3oo) < 4 MeV In the factorization 
approximation for the weak vertex we then get 

BR{Dt ^ P°7r+),(i3oo) < 7 x 10"^ , (3) 
where we have used fo^ = 230 MeV, to- 
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Fi gurG 1 . Annihilation diagram of Dg — > ptt decay. 



gether with the conservative assumptions of 
Bi?(7r(1300) pn) ~ 100% and r(7r(13GG)) 
equal to its upper experimental bound of 600 
MeV. The interference with other annihilation 
contributions from intermediate tt and 7r(1800) 
states can somewhat change the above estimate 
(using PCAC, the contribution from n was found 
in to be negligible, while the contribution 
of 7r(1800) is difficult to estimate due to the 
lack of experimental data). In addition, also the 
FSI contributions (to be considered shortly) fall 
in exactly the same range JHI- Therefore, un- 
less there are large cancellations, the value of 
BR{Df P^T^^) is expected to be near to its 
present experimental upper bound Q). 

In the case of the [ujtt^) final state there is 
no such resonance annihilation contribution and 
one has to explain a relatively large experimental 
value for BR{Ds lott^) in a different way. 
An important observation is that there are multi- 
body intermediate states that do have the correct 
values of J*^ and , for instance the two-body 
Ki*)K^*) states. As we will show in the rest of 
the talk, it is possible to explain the experimen- 
tal value for BR{Ds ojtt^) by considering the 
contributions due to the rescattering of these in- 
termediate states. 

In estimating the contributions from hidden 
strangeness intermediate states (that can arise 
from spectator quark diagrams), we use the fol- 
lowing assumptions 

• We consider only contributions coming 
from two body intermediate states with 
s, s quantum numbers (lowest lying pseu- 
doscalar and vector states). Note that the 
re-scattering through intermediate K, K* 



states is possible for both p'^7r+ as well as 
WTT"*" final state, while the re-scattering with 
intermediate rj or -q' is possible only in the 
case of WTT"*" final state due to isospin and G 
parity conservation. 

• For the weak transition D+ (X^*)/? (*') + 
in the D+ (/s:(*)x(*))+ p°Tr+ and 
D+ {K^*^K(*^)+ -> UJTT+ decay chains 
as well as for the weak transition Z?+ 
ri(ri')p^ in the Df — s- ri{ri')p'^ — > lutt^ decay 
chain we will use the factorization approx- 
imation. The weak Lagrangian is therefore 



X {ai{ud)H{sc)H + a2{sd)H{uc)H) , 

(4) 

with {ud)H,--- the hadronized V-A weak 
currents, Vij the CKM matrix elements and 
ai,2 the effective (phenomenological) Wil- 
son coefficients taken to be ai = 1.26 and 

02 = -0.52 mm\ . 

• Finally, the strong interactions are taken 
into account through the following effective 
Lagrangian |14I15I16] : 

(5) 

_4^^YR,f^-^0Tr{d^p,d^P0n) , 

where 11 and p^ are 3x3 matrices con- 
taining pseudoscalar and vector meson op- 
erators respectively and / is a pseudoscalar 
decay constant. We used numerical values 
Own = 0.33, and gprr. = 5.9 11 411 511 6| . 

In addition wc have checked that the use of fac- 
torization for the L»+ ^ K*+K*°, D+ K+K*° 
and Df K^K*^ decays gives reasonable es- 
timates of the measured rates (note that we do 
not need K'^K^ in further considerations) 

|l()j . In these results the annihilation contribu- 
tions have been neglected since they are an order 
of magnitude smaller. 
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r+ and D'^ —^ p'^7r+ decay amplitudes. 
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Figure 3. The intermediate rj, rj' , p'^ contributions ir 
the D'^ u)Tx^ decay. 
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Table 1 

Form factors at = [181 . The results in the first 
five columns are for Ds — » K,K*li/i transitions. The 
last column stands for the form factor appearing in the 
Ds — > rjslvi, (the ss component of rj, 77') transition. 



The situation in the case of 77, 77' intermediate 
states is not so favorable. To treat the 77, 77' mix- 
ing we use the approach of Ref. 19 with the value 
of the mixing angle transforming between 77, 77' 
and rjq ~ (mm + (iJ)/\/2, ?7s ~ ss states taken 
tohe <j> = 40°. The factorization approach then 
gives a reasonable description of £)+ p'^rj de- 
cay, while it does not reproduce satisfactorily the 
experimental result for —f p^Vf |10) . This is 
a known problem as the — > p'^r]' rate is very 
difficult to reproduce in any of the present ap- 
proaches |2l3l5j . This inevitably introduces some 
further uncertainty into our approach, yet the re- 
sulting uncertainty is not expected to affect sig- 
nificantly our main conclusions. 

For the weak current matrix elements between 
Ds and vector or pseudoscalar final states we 
use a common form factor decomposition I7ll()j 
with the form factors F(g^), Ai^2(q^)and 

Ao((?^). For the dependence of the form fac- 
tors we use results of ^S], based on a quark model 
calculation combined with a fit to lattice and ex- 
perimental data. Ref.^S] provides a simple fit 
to their numerical results with the form factors 
V{q^) and A^{c^') described by double 



pole dependence 



/(O) 



(l-q7M2)(l-crg2/M2) 
while single pole parameterization 
/(O) 



(1-CT(?7M2) 



(6) 



(7) 



can be used for ^1^2(9^), as the contributing res- 
onances have masses farther away from the phys- 
ical region (note that this parameterization ap- 
plies also to Fq form factor, which however does 
not contribute in the processes we discuss in this 
paper). The values of /(O) and a are listed in 
TablCnand are taken from [TH|. We use M = 1.97 
GeV in the expression for Aq, and M ~ 2.11 GeV 
for all the other form factors JS|. Incidentally, 
the parameterizations of the form factors © and 
(Q make all the loop diagrams in Figs. |21and|31 
finite. 

For the decay constants, defined through 
(0|q7''75(?|P(p)) = i!pv^ and (0|g7^g|F(p)) = 
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Table 2 

The dispersive Atn and absorptive AiA parts of the am- 
phtudes (in units of 10"'' GeV) for the — > wk^ decay 
corresponding to the diagrams on Fig. 13 (-Ai) and Fig. U] 
(Bn^n')- The amplitudes for the p^Tr"*" decay (ne- 

glecting the mass difference between rrxp and rriu) are ob- 
tained by inverting the sign of Ai£),AiA for even i, while 
= 0. 



gve^", we use fo = 0.207 GeV and fos = 1-13/r, 
as obtained on the lattice [TS] and for the rest 
fK = 0.16 GeV, \gK*\ = 0.19 GeV^ \gp\ = 
0.17 GeV^ and \g^\ = 0.15 GeV^ coming from 
the experimental measurements (71. 

The amplitudes for the _D+ con^ and _D+ 
p^TT"*" decays can be written as: 

A{Dt^c.n+)^^e-k,[j2^'r'+B), 

(8) 

A{Dt^p\+) = ^e-k,Y.4^\ (9) 

with e the hclicity zero polarization vector of the 
Lu or p vector mesons, while ^2 is the pion mo- 
mentum. The reduced amplitudes aI'''*''"'' and B 
correspond to the diagrams in Figs. [3 and 13 re- 
spectively. The explicit expressions can be found 
in Appendix of |10j . The numerical values for 
A[''^^''^^ and B are given in Table |21 Combining 
the above results we arrive at the prediction 

BR{Dt a>7r+) = 3.0 x 10"^ (10) 

Note that in this calculation we have used the 
factorization approximation for the diagram of 
Fig. 13 which as stated above, does not work 



well for Z)+ — > p+r/,?/ transition. Including 
hidden strangeness FSI to the Df p^rj' de- 
cay mode gives an order of magnitude smaller 
contribution. On the other hand one can use 
the experimental input to rescale the correspond- 
ing amplitudes. This results in the prediction 
BR{D+ W7r+) = 4.4 x IQ-^. We point out 
that the loop contributions are finite due to dou- 
ble pole parametrization of the form factors. If a 
single pole parametrization is used, one has to 
regularize the amplitudes. We found that the 
numerical results do not change significantly in 
this case when the cut-off scale is above but close 
enough to the Ds meson mass. We can draw 
the conclusion that the experimental result for 
BR{D'^ UJTT^) can be understood as a result 
of the combined effect of a spectator transition 
and FSI. Therefore, it makes the attempt to un- 
derstand the — > W7r+ amplitude as a result of 
annihilation contributions unsuccessful. 

In the case of the Df p^ir^ transition, the 
FSI contributions alone result in 

Si?(L>+ ^ p°7r+)FSi = 0.7 X 10"^ (11) 

This is almost exactly the same as our estimate of 
the upper bound on the annihilation contribution 
Q . Both contributions are equal or very close to 
the present 90% CL upper bound. If there is no 
destructive interference between these two con- 
tributions and the contributions of FSI through 
higher resonances that we did not take into ac- 
count, one hopes that the branching fraction for 
this decay will be determined in the near future. 
Our prediction is in agreement with the results of 
other theoretical studies which give the rate for 
£)+ p^n^ to be equal [3] or even larger than 
the rate for uir^ decay |5l6j . 

However, one should consider possible cancella- 
tion that might occur. Adding the FSI contribu- 
tion and the maximal annihilation contributions 
(PI with alternating signs gives a fairly large in- 
terval 

BR{Dt p"tt+) = (0.05 - 3.5) x 10"^ . (12) 

We note that the experimental uncertainties re- 
flected in the input parameters can change the 
values for BR{D+ ^ p°n+) and BR{D+ 
cj7r+) by about 20%. 



Finally, we mention that the kind of FSI con- 
tributions we were considering in this paper is 
not the leading contribution in the — > (/)7r+ 
transition, which can proceed through spectator 
quark transition directly. Use of the factorization 
approximation for the weak vertex leads to a pre- 
diction BR{D+ (/)7r+) = 4.0%, which is already 
in excellent agreement with the experimental re- 
sult of 3.6 ± 0.9%. We found that inclusion of 
FSI reduces the theoretical prediction from 4% to 

3.6%. The size of the shift also indicates that 
FSI of the type described in the present paper 
are in the case of Df — » 07r+ transition a second 
order effect. Note as well, that the size of the 
FSI correction is in agreement with the predic- 
tions for BR{D+ p°7r+) and BR{D+ a;7r+), 
which are an order of magnitude smaller than 
BR{D+ (j)TT+). 

We summarize that the hidden strangeness fi- 
nal state interactions are very important in un- 
derstanding the Df ujTT^ and Df p"7r+ de- 
cay mechanism. The Df — » wtt^ amplitude can 
be explained fully by this mechanism. In the case 
of the Df p^TT^ decay rate we obtain a fairly 
large range due to possible cancellation between 
FSI and single pole contributions. 

The measurement of the; Df ff'Ti'^ decay 
rate will considerably improve our understanding 
of the Df p^TT^ decay mechanism. The hidden 
strangeness FSI might fully explain the observed 
decay rate for BR{Df uin'^). Finally, this kind 
of FSI gives only subdominant contributions in 
the case of Dg (fnr, KK* decays, which are well 
described by the factorization approximation. 
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